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Abstract 

A  mathematical  modeling  of  a  system  consisting  of  a  cascade  of  a  thermionic  energy  conversion  (TIEC)  device  and  an  alkali  metal  thermal  to 
electric  converter  (AMTEC)  device  has  been  performed  to  optimize  the  efficiency  of  the  cell.  The  TIEC  is  heated  by  electronic  bombardment, 
which  converts  heat  partially  into  electricity  and  rejects  the  remaining  heat.  The  AMTEC  utilizes  the  rejected  heat  of  the  TIEC.  Cascading 
these  two  cells  provides  lots  of  advantages.  A  mathematical  model  for  the  cascade  converter  has  been  developed  to  analyze  the  effects  of 
key  parameters  such  as  power  level,  heat  fluxes  and  temperatures.  In  this  effort,  a  12-node  system  of  non-linear  simultaneous  equations  has 
been  constructed  which  is  solved  by  MATCAD  and  a  locally  optimized  efficiency  has  been  derived.  Thus,  efficiency  of  the  cascaded  cell  is 
improved,  so  that  it  is  greater  than  the  highest  efficiency  among  the  TIEC  and  AMTEC  and  lower  than  the  sum  of  their  individual  efficiencies. 
©  2005  Published  by  Elsevier  B.V. 
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1.  Introduction 

The  conversion  of  thermal  and/or  mechanical  energy  to 
electrical  energy  has  been  historically  the  mainstay  for  power 
systems.  With  the  advent  of  Faraday’s  law  of  electromagnetic 
induction  and  steam  engine,  dynamic  conversion  systems 
evolved  rapidly  and  were  perfected  with  the  engineering 
details  mostly  for  terrestrial  and  aircraft  devices.  As  a  poten¬ 
tial  alternative  to  the  dynamic  conversion  systems,  static  ther¬ 
mal  to  electrical  conversion  systems  have  been  investigated, 
particularly  for  space  power  systems.  With  the  discovery  of 
the  Seeback  effect,  the  thermo-electromotive  force,  which 
occurs  in  materials  under  the  effect  of  a  temperature  gradient, 
has  been  the  driving  force  for  static  thermo-electric  genera¬ 
tors.  A  large  number  of  static  devices  have  been  evolved  and 
studied  [1-9].  The  goal  for  static  conversion  systems  has  been 
to  identify  static  converters  with  efficiencies  that  are  compet¬ 
itive  with  dynamic  systems.  Two  compatible  static  devices 
may  be  linked  together  in  a  tandem  form  to  increase  the  over¬ 
all  efficiency.  The  objective  of  this  work  is  to  optimize  the 
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efficiency  of  a  device  that  cascades  two  static  devices,  namely 
TIEC  and  AMTEC. 

A  thermionic  energy  converter  (TIEC)  is  a  static  device 
(see  Fig.  1),  which  converts  heat  energy  through  the  surface 
emission  of  electrons.  It  consists  of  an  emitter  (cathode), 
which  receives  heat  from  a  suitable  source  and  emits  elec¬ 
trons  from  the  other  side.  A  collector  (anode)  collects  these 
electrons  and  is  cooled  to  a  lower  temperature  than  that  of  the 
emitter  to  limit  the  back  emission  of  electrons  [2].  The  anode 
and  the  cathode  are  connected  by  external  electrical  leads  to 
supply  the  generated  power  to  the  load. 

An  alkali  metal  thermal  to  electrical  converter  (AMTEC) 
is  a  device  for  the  direct  conversion  of  heat  to  electrical 
energy  with  no  moving  parts  (see  Fig.  2).  It  has  the  poten¬ 
tial  to  produce  electrical  power  with  much  higher  efficiency 
and  power  density  than  current  devices  such  as  the  radioiso¬ 
tope  thermo-electric  generator  (RTG),  general  purpose  heat 
source  modules  (GPHS)  [3].  AMTEC  uses  a  beta- alumina 
solid  electrolyte  (BASE)  to  separate  electrons  from  ions  in  a 
high-pressure  and  high-temperature  region  of  liquid  metal 
at  900-1300 K  [10].  It  converts  the  work  of  the  isother¬ 
mal  expansion  of  sodium  vapor  directly  to  electrical  power. 
AMTEC  has  many  advantages  for  terrestrial  and  space  power 
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Emitter, W 


Emitter  sleeve,  Mo/Re 


Collectror,  Ni 


The  input  for  TIEC  is  318  watt  and  the  efficiency  of  the 
cell  is  determined  by  the  equation: 

(°*Te-c).e 
n  = 

=  9.61% 

<n*Te.b)‘Aee+ea(Te4-Te4) 

where 

and 

Qeiectric "  (0*Te-c).e=5933  watt 
Qeiectric  “  (0*Te-b)*Aee=23.47  watt 
Qrade  "  eG(Te4-Te4)  =  38.206  watt 

Fig.  1.  Schematic  of  thermal  model  of  thermo-ionic  converter  (TIEC). 


AMTEC  heat  pipe 


Fig.  2.  Schematic  of  thermal  model  of  alkali  metal  thermal  to  electrical 
converter  (AMTEC). 


applications  including  no  moving  parts  with  the  resulting 
potential  for  low  maintenance  and  high  durability  and  effi¬ 
ciency. 

An  assembly  of  TIEC  and  AMTEC  has  been  cascaded 
with  a  transition  piece  between  the  collector  of  TIEC  and  the 
heat  input  region  of  AMTEC  (see  Fig.  3).  The  transition  piece 
is  surrounded  by  a  trim  collar  (heater/cooler),  which  can  be 
heated  by  a  resistance  heating  element,  or  cooled  by  water 
or  by  air,  in  order  to  maintain  the  desired  temperature  in  the 
transition  piece.  The  heat  rejected  by  the  collector  of  the  TIEC 
is  picked  up  by  the  heat  input  region  of  the  AMTEC  through 
the  transition  piece  and  thus  both  the  TIEC  and  AMTEC  are 
operated. 

Lodhi  et  al.  [1]  had  developed  a  model  by  writing  nine 
nodal  equations  (using  conservation  of  energy)  at  strategic 
locations  within  the  cascaded  cell.  The  number  of  nodal 
points  has  been  increased  from  9  to  12.  The  new  nodal  points 
included  were  at  the  insulator  between  the  transition  piece  and 
the  AMTEC  evaporator,  base  tube,  vapor  zone  (in  the  low- 
pressure  region  between  the  base  tube  and  the  condenser). 
The  variable  identified  in  this  cell  is  temperature  at  vari¬ 
ous  points  of  the  cell.  Mathematical  equations  for  the  energy 
balance  and  heat  fluxes  have  been  developed.  Through  this 
mathematical  model  the  optimized  efficiency  of  the  cell  is  to 
be  achieved.  When  we  are  cascading  two  cells,  the  efficiency 
of  the  cascaded  cell  is  expected  to  be  between  the  highest 
efficiency  of  the  individual  cells  and  less  than  the  sum  of  the 
efficiencies  of  the  two  cells.  For  this  process,  a  mathematical 
model  has  been  written  in  Mathcad,  Version  7.0. 


2.  Thermionic  energy  converter 

A  thermionic  energy  converter  is  a  device,  which  converts 
the  thermal  energy  to  electrical  energy.  There  are  two  types 
of  TIEC,  one  which  works  in  vacuum  and  the  other  works 
with  a  plasma. 

The  working  principle  of  vacuum  TIEC  involves  in  con¬ 
version  of  a  gas  of  free  electrons  in  the  inter-electrode  space, 
in  a  vacuum.  These  electrons  create  a  negative  space  charge 
between  the  electrodes,  and  as  a  result,  a  retarding  field  devel¬ 
ops  which  deflects  some  of  the  emitted  electrons  back  to  the 
cathode.  The  gas  chosen  in  TIEC  is  cesium  (Cs).  Cesium  has 
the  lower  ionization-potential  among  candidate  substances 
and  its  adsorption  reduces  the  surface  work  function  of  the 
electrodes  below  that  of  liquid  cesium  or  any  metal  [11]. 

The  cell  operates  at  a  temperature  of  1 800-2000  K  at  the 
emitter  end,  T  (see  Fig.  1).  It  has  been  designed  to  work 
independent  of  the  heating  source.  The  heat  from  emitter  zone 
is  either  radiated  or  conducted  to  the  other  parts  of  the  cell. 
A  portion  of  the  heat  is  radiated  to  the  surroundings,  Tam b. 
The  heat  gained  by  the  flange  from  the  emitter  is  in  two-fold, 
one  through  radiation  and  the  other  through  conduction.  The 
temperature  at  flange  is  7f.  The  gas  uses  some  amount  of  heat 
to  get  ionized  and  to  dissipate  the  electrons.  The  electrons 
get  accelerated  and  travel  through  the  electrodes  towards  the 
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The  input  for  TIEC  is  318  watt  and  the  efficiency 
of  the  cell  is  determined  by  the  equation: 


Qtrim  =  4.514  watt 


=  14.398% 


^electron  +^rade‘^trim 
Qelectric  =  (n‘Te-b)‘  Aee  =  23  47  watt 
Qelectric  =  (0*Vc)  e  =  5  933  watt 

Qrade  =  «*(Te4  '  Tc4)  =  38.206  watt 
^amtecelec  —  9  Ease'®  —  4.514  watt 

Qtrim  =  eWd*r(Tb4-Ttrim4) 

=  -10.878  watt 


Fig.  3.  Schematic  of  thermal  model  of  cascaded  TIEC/AMTEC  converter. 


collector.  The  electricity  is  drawn  from  the  electrodes  at  this 
point.  The  heat  energy  is  also  radiated  to  the  collector. 

The  heat  energy  from  the  flange  is  conducted  to  the  bel¬ 
lows,  radiated  to  surroundings  and  to  the  collector.  The  tem¬ 
perature  at  the  bellows,  7bei,  is  considerably  less  than  the 
temperatures  at  other  nodes.  Most  of  the  heat  is  given  out  to 
surroundings  from  bellows.  The  temperature  at  the  collector, 
Tc,  is  expected  to  be  between  1000  and  1250  K.  Release  of 
heat  energy  at  such  high  temperatures  makes  the  cell  ineffi¬ 
cient. 

The  converter  efficiency  of  TIEC  is  calculated  as  the  ratio 
of  the  contact  potential  difference  which  is  equal  to  the  volt¬ 
age  on  the  load  to  the  sum  of  the  energy  carried  away  by 
electrons  from  the  cathode  and  the  radiation  beat  transfer 
from  emitter  to  collector.  Thus,  the  efficiency,  ??tiec>  is: 


77TIEC  = 


Q  electric 


2  electron  T  2rade 


where  2eiectric  =  (0  x  Te  —  c)  x  Aee,  O  =  constant  =  0.0135 
WK-1cm-2,  Te  =  emitter  temperature  in  K,  Aee  =  area  of 
the  emitter  in  cm2,  c- power  per  unit  area  in  Won-2; 
2electron  =  (nxTe  —  b)  x  Aee,  n  =  constant  =  0.0423  W  K~ 1 
cm-2,  £>  =  58.6 Won-2;  2rade  =  <J  x  s  x  (Te4  —  Tc4), 

<7  =  5.61  x  10-12  Won-2  K~4,  e  =  emissivity  =  0.18,  Tc  = 
collector  temperature  in  K. 

3.  AMTEC 

AMTEC  is  a  high-temperature  regenerative  concentration 
cell  for  an  appropriate  alkali  metal,  which  converts  thermal 
energy  directly  into  electrical  energy  (see  Fig.  2).  In  this 
study,  liquid  sodium  is  selected  as  the  working  fluid.  The 
efficient  operation  of  the  AMTEC  cell  involves  several  heat 
and  mass  transfer  processes.  AMTEC  is  a  relatively  new 
type  of  device,  based  upon  the  principle  of  an  alkali  metal 
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concentration  cell,  conceived  in  late  1960s.  This  is  a  static 
conversion  device,  which  can  provide  efficiency  close  to  the 
theoretical  Carnot  efficiency.  AMTFC  operates  at  a  temper¬ 
ature  of  900-1300  K  at  hot  side  and  400-600  K  at  the  cold 
side. 

There  are  two  types  of  AMTEC  cells  based  on  their  work¬ 
ing  principle.  One  is  liquid  anode  type  and  the  other  is  vapor 
anode  type.  In  this  study,  we  are  simulating  the  liquid  anode 
type.  To  understand  the  working  of  AMTEC,  its  operating 
cycle  is  illustrated  in  Fig.  2  schematically.  A  closed  ves¬ 
sel  is  divided  into  a  high-temperature  (Ta)  region  in  contact 
with  a  heat  source  and  a  low-temperature  (7j)  at  the  cooler 
region.  The  critical  material  in  the  operation  of  the  AMTEC 
is  the  beta-alumina  solid  electrolyte,  a  sodium  ion  conductor 
whose  ionic  conductivity  is  much  larger  than  its  electronic 
conductivity.  In  an  AMTEC  cell,  the  BASE  separates  the 
hot  (high-pressure)  region  filled  with  a  small  quantity  of  liq¬ 
uid  sodium  in  contact  with  the  heat  source,  from  the  cold 
(low-pressure)  region  occupied  with  sodium  vapor.  A  porous 
electrode  covers  the  low-pressure  (outer)  side  of  the  BASE. 
Electrical  leads  in  contact  with  the  porous  electrode  and  the 
high-temperature  liquid  sodium  exit  through  the  wall  of  the 
device  and  is  connected  to  an  external  load.  The  pressure  dif¬ 
ferential  across  the  BASE  forces  ionization  of  sodium  atoms 
on  the  hot  side.  The  ions  diffuse  through  the  BASE  towards 
the  low-pressure  side  in  response  to  the  pressure  differential 
gradient  of  free  Gibbs  energy  while  the  electrons  circulate 
through  the  external  load,  producing  electrical  work.  Elec¬ 
trons  and  sodium  ions  recombine  at  the  interface  between  the 
BASE  and  the  porous  electrode.  The  resulting  sodium  ions 
absorb  their  heat  of  vaporization,  move  through  the  electrode 
and  the  vapor  space,  then  release  their  heat  of  vaporization 
at  the  low- temperature  side  of  the  condenser  surface.  A  wick 
structure  or  an  electromagnetic  pump  brings  the  liquid  con¬ 
densate  to  the  high-temperature  side  of  the  BASE  tube  to 
complete  its  circulation. 

The  heat  distribution  in  the  AMTFC  is  as  follows:  the 
temperature  at  the  hotplate,  ra,  should  be  between  900 
and  1300  K.  The  heat  source  for  AMTEC  may  be  radioac¬ 
tive  materials,  solar,  heat  rejected  from  other  devices  or 
any  other  heating  source.  The  heat  energy  from  this  point 
is  given  out  to  the  sodium  liquid.  Because  of  convective 
heat  transfer,  the  sodium  starts  boiling  and  slowly  gets  con¬ 
verted  to  vapor.  The  heat  energy  allows  the  sodium  to  move 
towards  the  low-pressured  BASE.  As  the  sodium  approaches 
the  BASE  it  gets  ionized.  The  temperature  at  the  BASE, 
7base,  is  approximately  equal  to  that  of  the  boiling  point 
temperature  of  sodium.  Hence  some  part  of  heat  energy  is 
used  in  transporting  sodium,  some  in  generating  electricity 
and  some  of  the  heat  energy  is  wasted  as  parasitic  losses. 
Also  some  part  of  heat  energy  is  transferred  to  the  sodium 
vapor  in  convection  mode.  The  heat  passed  onto  the  vapor 
zone  from  the  BASE  is  linearly  distributed  to  the  condenser 
and  to  the  surroundings  through  the  wall  of  the  cell.  The 
temperature  at  the  condenser,  7),  is  expected  to  be  above 
400  K. 


The  efficiency  of  AMTEC,  rj amtec  ,  is  defined  as  the  ratio 
of  the  amount  of  electric  energy  generated  to  the  amount  of 
heat  supplied  to  the  AMTEC. 

Thus,  the  efficiency  of  AMTEC  is: 

Q  amtecelec 

7?  AMTEC  =  — — - 

Grade 

where  gamtecelec  =  q  x  7base  -  e,  q  =  constant  =  1.2  x  10“ 12 
WK-1,  7|laSL.  =  BASE  tube  temperature  in  K,  e  =  8.9 W; 
grade  =  a  x  e  x  (if  -  r4),  or  =  5.67  x  10-12  W  cm-2  K-4, 
s  =  emissivity  =  0.18,  Tc  =  collector  temperature  in  K. 

4.  Cascade  of  TIEC  and  AMTEC  cell 

The  purpose  of  cascading  TIEC  and  AMTEC  is  to  derive 
an  optimum  power  output  with  less  heat  energy  wasted.  A 
schematic  diagram  of  such  a  cascade  is  shown  in  Fig.  3.  Cas¬ 
cading  of  direct  energy  conversion  devices  has  been  designed 
for  the  waste  heat  of  a  high-temperature  device  used  to  oper¬ 
ate  a  bottoming  low-temperature  device.  It  would  allow  the 
development  of  a  highly  efficient,  compact  and  light  weight 
power  to  address  future  civilian  and  defense  missions  [5]. 
This  system  uses  a  cesium-TIEC  to  operate  a  bottoming 
sodium- AMTEC  [12].  This  combination  is  most  suitable 
because  the  optimum  collector  temperature,  Tc ,  of  TIEC  is 
typically  in  the  1000-1200  K  ranges  and  AMTEC  hot  pipe 
temperature,  ra,  falls  well  within  this  range. 

The  construction  of  such  system  involves  using  a  transition 
piece  between  the  collector  of  TIEC  and  the  heat  pipe  of 
AMTEC.  The  transition  piece  is  surrounded  by  a  trim  collar 
(heater),  which  can  be  heated  by  a  resistance  heating  element, 
or  cooled  by  water  or  air,  to  maintain  the  desired  temperature 
and  heat  flux  in  the  transition  region.  The  TIEC  is  heated  up 
by  the  electron  bombardment  on  the  upper  end  of  the  emitter, 
which  converts  a  portion  of  heat  energy  into  electricity  and 
rejects  the  rest  from  the  lower  end  of  the  collector,  into  the 
transition  piece.  The  remaining  heat  rejected  at  the  bottom  of 
the  collector  is  picked  up  by  the  AMTEC  via  the  transition 
piece. 

The  heat  rejected  from  the  collector  of  TIEC  is  absorbed 
by  the  transition  piece.  Some  portion  of  the  heat  energy  may 
be  used  to  heat  the  trim  collar.  The  heat  energy  supplied  to 
trim  collar  will  be  through  radiation.  If  an  excess  amount 
of  heat  is  supplied  to  the  transition  piece,  then  some  heat 
is  passed  onto  the  trim  collar.  If  the  collector  of  TIEC  does 
not  supply  sufficient  amount  of  heat  energy  to  the  transition 
piece,  the  trim  collar  will  supply  the  required  heat  energy. 
Hence,  the  heat  contribution  by  the  trim  collar  to  the  total 
cell  may  be  included  or  excluded.  A  calorimeter  is  used  for 
the  transition  piece  to  measure  the  amount  of  heat  energy 
passing  through  it.  Thus,  certain  amount  of  heat  energy  is 
radiated  from  the  trim  collar  to  the  calorimeter,  and  part  of  it 
from  the  calorimeter  to  the  surroundings.  A  certain  amount 
of  heat  energy  is  supplied  to  the  AMTEC  heat  pipe.  The  heat 
flow  of  the  cascaded  cell  is  illustrated  in  Fig.  4. 
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Fig.  4.  Thermal  circuit  for  calculating  temperatures  and  power  flows  in  thermo-ionic/ AMTEC  cascade  converter. 


The  efficiency  of  the  cascaded  cell  should  be  larger  than 
the  highest  efficiency  of  the  individual  cells  and  lesser  than 
the  sum  of  the  efficiencies  of  both  the  cells.  The  efficiency  of 
the  cascaded  cell,  r] cascade*  is: 


^cascade  — 


2  electric  T  2  amtecelec 
2  electron  T  2rade  —  2  trim 


where  2electric  =  (0  xTq  —  c)  x  Aee,  O  =  constant  =  0.0135 
WK-1cm-2,  re  =  emitter  temperature  in  K,  Aee  =  area  of 
the  emitter  in  cm2,  c-  power  per  unit  area  in  Wcm-2; 
Qamtecelec  =  <7  X  7baSe  “  e,  Cj  =  constant  =  1.2  X  10-12  WK"1, 
^base  —  BASE  tube  temperature  in  K,  e  =  8.9W;  2eiectron  = 
(n  x  Tq  —  b)  x  Aee,  n  =  constant  =  0.0423  WK-1  cm-2,  b- 
58.6 Wcm-2;  2rade  =  cr  x  s  x  (T4  —  T4),  a  =  5.61  x 
10-12  Wcm-2  K-4,  £  =  emissivity  =  0.18,  Tc  =  collector 

temperature  in  K;  2trim  =  cr  x  s  x  n  x  <7trans  x  /trans  x 
(T^  —  r4im),or  =  5.67  x  10-12  Wcm-2  K-4,  £  =  emissivity  = 
0.2,  dtrans  =  diameter  of  the  transition  piece  in  cm, 
^trans  =  length  of  the  transition  piece  in  cm,  rtrim  =  trim 
collar  temperature  in  K,  T t,  =  transition  temperature  in  K. 

The  mathematical  model  is  written  for  the  cascaded  cell 
and  is  aimed  at  optimizing  the  efficiency  of  the  cell. 


5.  Mathematical  model  and  feasible  solution 

A  mathematical  model  has  been  developed  in  order  to 
maximize  the  efficiency  of  the  cascaded  cell  satisfying  all  the 
constraints  associated  with.  The  cascade  cell  has  been  divided 
into  12  nodal  points  and  at  each  nodal  point  the  heat  balance 


equations  have  been  written.  These  heat  balance  equations  are 
the  constraints  to  the  problem.  As  a  general  rule  heat  always 
travel  from  hot  region  to  cold,  hence  the  other  constraints  of 
the  problem  are  Te  >  Tf  >  Tc  >  Tb  >  Ta  >  7base  >  Tsodv  >  T\. 

The  objective  function  is: 


Maximize 


2  electric  T  2  amtecelec 
2electrn  T  2rade  —  2  trim 


The  heat  balance  constraints  are  given  in  Appendix  A.  The 
constraints  listed  in  Appendix  A  are  the  12  nodal  points  which 
are  located  at  the  emitter,  flange,  bellows,  collector,  transition 
piece,  trim  collar,  calorimeter,  insulator,  AMTEC  evapora¬ 
tor  (at  the  high-temperature  region),  base  tube,  vapor  zone 
(low-pressure  region)  and  condenser  (at  the  low-temperature 
region),  respectively.  At  each  of  these  points  a  detailed  study 
of  heat  distribution  is  done  and  the  resulting  equations  are 
constrained  to  the  problem. 

Apart  from  these  constraints,  there  are  other  parameters, 
which  significantly  affect  the  efficiency  of  the  cell.  These 
terms  are  involved  in  AMTEC  for  parasitic  losses,  electric 
power  output  and  sodium  heat  transportation.  These  parame¬ 
ters  affect  the  AMTEC  cells  efficiency,  and  thus  the  efficiency 
of  the  cascaded  cell.  Once  the  program  has  been  run  and 
checked  for  feasible  solution,  these  variables  are  included 
and  again  the  program  is  run.  Thus,  by  applying  the  search 
method  for  the  problem  with  these  constraints  a  local  optimal 
solution  is  achieved. 

The  objective  function  and  the  constraints  are  pro¬ 
grammed  in  Mathcad  and  are  allowed  to  run  for  several 
iterations  with  a  minimum  error  factor  of  1.507  x  10-12. 
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Table  1 


Temperature  profile 


Norder 

Temperature  (K) 

Te 

1940 

A 

1310 

Tbel 

1036 

Tc 

1225 

Tb 

1184 

Ttrim 

1229 

Tins 

1167 

Ta 

1159 

Tbase 

1118 

Tsodv 

635.941 

Ti 

413.746 

Teal 

878.045 

The  constraints  for  the  temperature  profile  are:  (1)  Te>Tf>Tc>T\j> 
T&  >  7base  >  Tsodv  >T\\  (2)  1100K>  Tbase  >  1 150 K;  (3)  400 K<7i  < 500 K; 
(4)rf>rbei;(5)re<2000K. 


7.  Conclusion 

The  basic  purpose  of  optimizing  the  cascaded  cell  is 
to  limit  weight,  area  and  material  property  constraints  and 
finally  improve  the  cost  effectiveness.  These  characteristics 
have  been  achieved  by  predicting  a  temperature  profile  in 
Table  1 .  The  efficiency  of  the  cascaded  cell  has  been  obtained 
as  14.398%.  Though  the  solution  is  locally  optimized  still  fur¬ 
ther  work  needs  to  be  done  to  globally  optimize  the  efficiency. 
This  project  is  a  good  beginning  towards  the  optimization  of 
the  cascaded  cell.  In  this  direction,  the  Mathcad  program  is 
an  acceptable  approach  for  optimization  of  the  cell  to  the 
extent  that  it  could  not  handle  any  more  equations  than  used 
in  this  program.  For  a  program  containing  more  equations 
than  the  present  work  does  need  some  other  approach  with 
larger  space  and  capability. 


Mathcad  program  simultaneously  solves  for  the  temperature 
profile  and  the  search  method.  Each  time  the  program  runs  in 
iterations  and  solves  the  problem,  the  aforementioned  param¬ 
eters  are  changed  in  an  order.  Thus,  a  feasible  solution  is 
achieved  which  satisfies  all  the  constraints.  The  model  had 
given  a  feasible  solution  of  14.398%  satisfying  all  the  con¬ 
straints.  This  solution  of  14.398%  is  the  efficiency  of  the 
cascade,  which  is  greater  than  the  highest  of  the  individ¬ 
ual  cells  where  the  efficiencies  of  individual  cells  are,  for 
TIEC  9.66%  and  for  AMTEC  11.81%  and  lesser  than  the 
sum  of  the  efficiencies  of  the  both  the  cells,  21.47%.  Hence, 
the  solution  is  acceptable.  The  temperature  profile  is  given  in 
Table  1. 


6.  Result  and  discussion 

The  Mathcad  program  gives  a  locally  feasible  solution 
of  the  mathematical  model  constructed  in  this  work.  For 
any  unconstrained  non-linear  program  there  can  be  several 
optimal  solutions,  and  for  the  maximization  problem,  among 
these  optimal  solutions  again  there  can  be  a  maximum  val¬ 
ued  solution.  This  maximum  valued  solution  is  known  as  the 
global  optimal  solution  and  the  rest  of  them  are  the  local 
optimal  solutions.  In  this  problem  also  a  solution  has  been 
achieved  and  this  solution  may  or  may  not  be  the  global  opti¬ 
mal  solution.  It  is  found  that  the  efficiency  of  the  TIEC  is 
9.66%  and  that  of  the  AMTEC  is  11.81%.  The  efficiency 
of  the  cascaded  cell  is  obtained  to  be  14.398.  These  results 
have  satisfied  all  the  constraints  defined  in  the  mathematical 
model.  The  heat  balance  equations  have  been  formed  and  are 
allowed  to  run  in  iterations,  and  the  result  showed  that  the 
temperature  profile  is  satisfactory.  The  temperature  profile  is 
given  in  Table  1 . 

The  results  achieved  by  this  project  were  comparatively 
better  than  the  results  of  the  previous  work,  which  included 
nine  nodal  points  [1]. 
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Appendix  A 


For  the  objective  function, 

A  j  Q  electric  T  Q  amtecelec 

Maximize 


Gelectrn  T"  Grade  Gtrim 
The  heat  balance  constraints  are: 


1. 

2. 

3. 


5. 

6. 

7. 

8. 

9. 

10. 

11. 

12. 


16.13  x  10-12  x  T4  -  9.348  x  10“12  x  r4mb  -  3.206  x  10~12  x 
T4  -  3.576  x  10-12  x  If  +  0.1792  x  Te  -  0.1234  x  Tf  <  239.14 
3.576  x  10-12  x  T4  -  32.022321  x  10“12  x  T4  +  6.5225208  x 
10-12  X  r4el  +  15.154  X  10-12  X  r4mb  +  6.7688  x  10“12  x  T4  + 
0.123  x  Te  -  0.22987f  +  0.01  x  Tamb  +  0.08  x  Tbei  +  5.73  <  0 
6.5225208  x  10“12  x  Tf  -  69.210247  x  10“12 


x  T4el  +  37.106  x 


10-12  x  T4mb  +  0.086  x  7f  +  0.086  x  Tc  -  0.172  x  Tbei  + 
4.153126  x  10-12  x  T4  +  21.4486  x  10“12  x  T4al  <  0 
4.  4.1626  x  10-12  x  T4  -  16.7148  x  10“12  x  T4  +  6.7688  x 
10-12  x  T4  +  5.7834  x  10~12  x  T4mb  +  0.0423  x  Te  +  0.086  x 


rbe i  -  0.9887  x  Tc  +  0.9072  x  Tb  <  58.6 

0.9027  x  Tc  -  3.1807  x  Th  +  2.278  x  71ns  -  43.673955  x  10“12  x 
Tf  +  34.476377  x  10“12  x  Tfm  +  9.197578  x  10“12  x  T4el  <  0 


34.476377  x  10' 


-12 


xTb4 

46.82815  x  1(T12  x  r4al  +  90  <  0 


trim 

81.304527  x  10“ 


■12 


x  Tf  + 

trim  1 


21.4486  x  10 
10 


i—12  w  t^4 


12  x  4, 


222.3387  x  10 
x  T 


-12  --4  <0 


12  x  rc4al  +46.8281  x 


x  7^  +  154.06196  x  10 

1.8224  x  Tb  -  7.2992  x  rins  +  34.476377  x  10“12  x  714 


trim 


37.556844  x  10“12  x  T.4  +  5.4772  x  ra  +  2.7930738  x  10“12  x 

ms 

Cmb  <  0 

5.4772  x  7ins  x  T2  +  5.4772  x  71ns  x  T2d  +  10.9544  x  7jns  x 
Tsod  X  Ta  -  291829.48  x  Ta3  +  875499.39  x  T2od  x  Ta  - 
875493.9  x  Tsod  x  T2  -  291829.48  x  T3od  <  0 
291829.48  x  Ta3  -  875488.44  x  T2  x  rsod  +  875488.44  x  Ta  x 
Cod  “  291829.48  x  T3od  -  2.068977  x  Tbase  +  0.0068977  x 
Tsodv  +  53.2  <  0 

0.006977  x  Tba se  -  0.0197813  x  rsodv  +  0.0058273  x  Tamb  + 
0.006977  x  7]  <  0 

2.05  x  Tba se  +  0.006977  x  Tsodv  -  0.406977  x  Tx  +  0.4  x  Tamb  < 
44.3 
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The  constraints  listed  above  are  the  12  nodal  points,  which 
are  located  on  the  cascade.  At  each  of  these  points  a  detailed 
study  of  heat  distribution  is  done  and  the  resulting  equations 
are  constrained  to  the  problem. 
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